Normal Helium 3: a Mott-Stoner liquid. 
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A physical picture of normal liquid 3 He, which accounts for both "almost localized" and "almost 
ferromagnetic" aspects, is proposed and confronted to experiments. 
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Helium 3 is a liquid of strongly interacting fermionic 
atoms. As pressure is increased between p = and 
34bars, the compressibility is drastically reduced from 
k/kq ~ .27 to k/kq ~ .066, while the effective mass (spe- 
cific heat coefficient) and magnetic susceptibility are en- 
hanced from m*/m ~ 2.8, x/xo — 9.2 to m*/m ~ 5.8, 
x/xo — 24 (m is the mass of free 3 He atoms, and kq, Xa 
refer to the free Fermi gas at the same pressure) Q . This 
is clearly the signal of the increasingly important effect 
of the interatomic interaction, which has a strong repul- 
sive hard-core. Given the difficulties of a full quantita- 
tive treatment of such a problem, more phenomenological 
theoretical descriptions have been sought. 

There are mainly two different, and seemingly contra- 
dictory, physical pictures that have emerged over the 
years. In the "almost ferromagnetic" approach 
the liquid is viewed as being increasingly close to a fer- 
romagnetic instability as pressure is increased and the 
large susceptibility is explained as a Stoner enhancement 
X = Xo/0- — IXo)- Being very close to an instability, crit- 
ical spin- fluctuation modes ( "paramagnons" (PM)) must 
be taken into account beyond the mean-field Stoner de- 
scription. These modes have been claimed to be essential 
to explain, for example, the low-temperature dependence 
of the susceptibility [Q . They also provide a logarithmic 
increase of the effective mass (though in rather mediocre 
quantitative agreement with experimental values). The 
"almost ferromagnetic" PM picture has some severe lim- 
itations however, particularly in failing to explain the 
strong reduction of compressibility || . 

This reduction is one of the main motivation for view- 
ing instead the liquid as "almost localized" , i. e becoming 
more and more "solid-like" with pressure, as first pro- 
posed by Anderson and Brinkman g), and extensively 
developed by Vollhardt || . The essential physics behind 
that picture is that of localization by repulsive interac- 
tions, in the sense of Mott. As explained below however, 
the simplest implementation of the quasi-localized pic- 
ture leads to an incorrect description of the magnetic 
correlations and spin-fluctuation properties. 

In this letter, we would like to propose a novel descrip- 
tion of liquid 3 He, which retains the proximity to Mott 



localization as a central notion, but reintroduces a more 
accurate description of the spin-spin correlations, of dom- 
inantly ferromagnetic nature. In our picture, in addition 
to being "almost localized", the liquid is also close to a 
ferromagnetic instability (see also 0), but not in a crit- 
ical regime (contrary to PM theory). For this reason, 
liquid 3 He is viewed in our picture as a "Mott-Stoner" 
liquid. 

The quasi-localized picture was first implemented 
quantitatively [^|J^] by considering a lattice-gas model 
and modeling the hard-core repulsion as a Hubbard in- 
teraction: 
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This model was then treated using the Gutzwiller ap- 
proximation (GA), which yields a Mott transition at half- 
filling [n = 1) between a Fermi- liquid phase for U < U c 
and a localized Mott insulating phase for U > U c . Close 
to the transition, the compressibility vanishes as U c — U 
and the effective mass diverges as 1/(U C — U), in qualita- 
tive similarity with the behavior of liquid 3 He for increas- 
ing pressure. In the simplest description, the lattice-gas 
is constrained to be at half-filling, and the GA expression 
for m*/m can be fitted to the experimental result in or- 
der to extract the single parameter u(p) = U(p)/U c , with 
u(p) an increasing function of p. (More refined formula- 
tions la] also allow a variable filling factor, and consider 
a trajectory u = u{p), n — n(p) in the [u, n) plane). 
Given u(p), the calculated GA compressibility is found 
to be in reasonable agreement with experiment. It is 
interesting in this respect to consider the dimensionless 
ratio R K = (k/ Ko){m* /m), which is predicted to reach 
a finite value at the transition R^ A (U C ) ~ .25. Exper- 
imentally, this ratio does saturate at high pressure at a 
value R e K xp ~ .38. 

Turning to magnetic properties, the GA also leads to a 
divergent susceptibility x ~ 1/(U C — U), and thus to a fi- 
nite "Wilson ratio" R w = 1/(1 + F a ) = {x/Xo)/[m* /m). 
Experimentally, this ratio has a very weak dependence on 
pressure, with R C w P {v = 34bars) ~ 4.1, close to R^ A (U C ). 
This agreement was originally viewed flq] as one of the 
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main success of the approach and interpreted as evidence 
that the susceptibility enhancement could be entirely due 
to the incipient localization S, responsible for the ef- 
fective mass enhancement. However, we would like to 
point out that this divergence of the uniform susceptibil- 
ity is an artefact of the GA, rather than a genuine fea- 
ture of the Hubbard model. Indeed, in this model, the 
superexchange mechanism produces a nearest-neighbor 
antiferromagnetic exchange (of order J ~ e 2 F /2U for 
large enough U). On physical grounds, one expects this 
magnetic exchange to cutoff the divergence of the uni- 
form susceptibility, which should remain finite, of order 
X — 1/J through the transition and in the localized 
phase. More accurate treatments of the Mott transi- 
tion within a dynamical mean- field of the Hubbard model 
based on the limit of large lattice coordination fully sup- 
port this view . Given the Fermi energy of liquid 3 Hc, 
the magnetic superexchange would induce antiferromag- 
netic correlations on the scale of J ~ 350mK, which is 
physically unrealistic and yields a too small susceptibil- 
ity enhancement x/xo — £f/J < 7. This has actually an 
even more drastic consequence, namely that the ground- 
state of the half- filled Hubbard model is in fact an antifer- 
romagnetic insulating solid, rather than a paramagnetic 
liquid (except if a very large lattice frustration is intro- 
duced JhJ). In the GA, the superexchange is neglected 
altogether so that these difficulties are simply overlooked. 
The overestimate of short-range antiferromagnetic corre- 
lations is actually not entirely due to the lattice descrip- 
tion of the system. Variational treatments in continuum 
space using Jastrow-Slater wave functions and a realistic 
interatomic potential, also suffer from similar problems 
(when confronted, e.g to neutron data) pjJTl]- This is 
because, as in the Hubbard model, the emphasis is put 
mainly on only one of the effects of the hard-core, namely 
the avoidance of double occupancy. 

Here, we suggest that the original formulation of the 
quasi-localized picture must be modified in order to ac- 
count for the correct scale and type of magnetic cor- 
relations in the liquid. A very simple way to achieve 
this, while remaining in the framework of a lattice-gas 
model as above, is to introduce explicitly an additional 
nearest-neighbor ferromagnetic exchange, leading to the 
two-parameter model: 

ij,a i <ij> 

The overall magnetic exchange 1=1— J~J — ej?/2U, 
can have a priori an arbitrary sign, and will be deter- 
mined below from experimental considerations. We shall 
use the dynamical mean-field approximation, formally 
exact in the limit of large lattice coordination z — > oo, 
to analyze the behavior of this model (for an extensive 
review and technical details see Ref. @). In this limit, 
the inter-site magnetic interaction can be treated at the 



(static) mean-field level (in contrast, the local Hubbard 
interaction requires a full dynamical treatment). Hence, 
the dynamical magnetic susceptibility x(<f, w) is readily 
obtained from that of the Hubbard model Xh{<1i at 
the same value of U as: 

Xiq^)- 1 = XHiq^r 1 - IA(q] (3) 

where A(q) is the Fourier transform of the 
nearest- neighbor connectivity matrix: A(g) = 

i Ei=i exp in- Rj)- 
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FIG. 1. Phase diagram of model (g). The squares indi- 
cate the fit to liquid 3 He data as pressure is increased from 
p = 0, 3, 6, • ■ • , 33bars. Also indicated by arrows are the cor- 
responding trajectories in the "almost localized" (GA) and 
the "almost ferromagnetic" (PM) pictures. 

In Fig.l, we display the zero-temperature phase di- 
agram of this model at half-filling, in the (£/, /) plane 
0,0. We shall mainly focus in the following on the 
Fermi liquid phase which is found when the Hubbard re- 
pulsion is smaller than a critical value associated with 
Mott localization (U < U c ~ 3.7^), and the ferromag- 
netic exchange / is in the intermediate range lf F (U) < 
I < I F (U). We emphasize that this phase does not dis- 
play any kind of magnetic long-range order (the overall 
magnetic exchange / being in an intermediate coupling 
regime, it is always successfully opposed by kinetic energy 
effects). Stabilizing such a phase is one of the primary 
motivation of our approach. (Note that, as expected, the 
half-filled pure Hubbard model (/ = 0) is ordered an- 
tiferromagnetically for all values of U at T = in our 
treatment). When / is too small (/ < I^ F (U)), the an- 
tiferromagnetic superexchange induced by the Hubbard 
repulsion takes over, and the liquid orders antiferromag- 
netically, while for large I > I F (U), it orders ferromag- 
netically. The Mott localized state for U > U c is always 
magnetically ordered at T = 0, either antiferromagneti- 
cally for / < 1^ or ferromagnetically for I > I F , with 
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I^ F (U) = if (U) = J ~ 4/217 for [/ > U c . 

Simple estimates of the critical couplings I F ,lf F can 
be obtained for small U, and U close to U c To first or- 
der in U, the Hubbard model static susceptibility reads: 
XH(q = 0)- 1 = Xo 1 -U = e F -U, and xniQ)- 1 = U 
for the antiferromagnetic wavevector Q = (tt, ■■■,tt). 
Inserting those into Eq.(||), one obtains for small U: 
If = e F -U + 0(U 2 ) and I^ F = U + 0{U 2 ). Close 
to the Mott transition, the static susceptibility of the 
Hubbard model in the dynamical mean-field approach is 
well approximated by the form p0| : xh(<1,u — 0) _1 ~ 
X(q)e F + JA(q). In this expression, X(q) has a rather 
weak ^-dependence (with A(0) = 1), while e* F is the low- 
energy effective Fermi scale e* F ~ .66Ze F , with Z the 
quasi-particle residue which vanishes as the Mott transi- 
tion is reached at half-filling: Z ~ .9(1 — U/U c ). Hence, 
the static (/-dependent susceptibility of (||) is reasonably 
approximated for U ~ U c by: 

X (q,U = 0)' 1 ~ X(q)e* F + (J - I) A(q) (4) 

Using this expression, one obtains the estimates for U 
close to U c : If ~ J + e* F ~ J + .16(C/ C - H), I c Ai " ~ 
J-A(Q)e>. 

We now focus on the behavior of various physical quan- 
tities in the strongly correlated liquid phase, and on the 
application to the physics of liquid 3 He. In this phase, 
the single-particle self-energy, and the density-density re- 
sponse function are unaffected by the coupling I within 
the dynamical mean- field treatment, so that the effective 
mass and compressibility are functions of U only and 
coincide with those of the pure Hubbard model. Close 
to the Mott transition, the effective mass diverges as 
m*/m= 1/Z~ 1.1/(1 — U/U c ), while the compressibility 
vanishes as k/kq — .66(1 — U/U c ), which is qualitatively 
similar to the GA (in contrast with \ — l/( e F + J ~ -0)- 
Following the original phenomenological spirit of the 
quasi-localized approach |6) , one can determine the pres- 
sure dependence of the two effective parameters U (p) , 
I(p) by fitting the experimental results for two physi- 
cal quantities. We have chosen to use the effective mass 
(specific heat) data to determine U(p), and then to ex- 
tract I(p) from the susceptibility. For each pressure, liq- 
uid 3 He thus corresponds to a point indicated on the 
phase diagram of Fig.l. The resulting trajectory is seen 
to approach both Mott localization U = U c and the fer- 
romagnetic phase boundary I — if (U) . While the 
effective mass enhancement is entirely associated with 
the on-site repulsion U, the proximity of the ferromag- 
netic phase boundary is crucial to account for the ob- 
served magnitude of the susceptibility enhancement: at 
the value of U corresponding to the highest pressures, the 
susceptibility enhancement of the pure Hubbard model 
would be xh/xo — 5.4, about 4.5 times too small. In 
the present approach, the saturation of the Wilson ratio 
Rw = 1/(1 + F§) precisely reflects the fact that both in- 



stabilities are approached. Indeed, using Eq. (0) we have, 
close to U c : x/xo = (f/{^ f + J - I) = e F /{lf - I), so 
that Rw = (x/xo)/{ m * l m ) — e F/(lf ~ I)j an d hence 
U c — U(p) ~ 4Rw{Ic — I{p)), indicating a linear tra- 
jectory towards the multicritical point at the highest 
pressure. The total magnetic exchange I = I — J ~ 
(.66Rw — l)/x — 1-6/x obtained from our approach 
is ferromagnetic (I > 0) and of the order of 300mK. 
This is precisely the typical energy gained by including 
short-range magnetic correlations in variational calcula- 
tions |l],[ll[]. The ferromagnetic sign is consistent with 
an instability of the liquid towards a triplet superfluid 
phase at low temperature. Having determined U(p), I{p), 
we have compared the compressibility computed for our 
model to experiment. Excellent agreement is found at 
low pressure, while calculated values at high pressure 
are too large by approximately a factor of 2. The ra- 
tio R K = (k/ 'Ko)(m* /m) is predicted to saturate at high 
pressure as observed experimentally (with R K (U C ) ~ .73, 
while R% xp (p = 34bars) ~ .38). 

We would now like to compare and contrast the phys- 
ical picture proposed here to that of the "almost ferro- 
magnetic" PM description We first evaluate the 
dimensionless parameter r = (if — I) /if measuring the 
distance to the ferromagnetic critical boundary. From 
above, we find r ~ \.b/Rw, which varies from r ~ .46 
at low pressure to r ~ .36 at high pressure and is thus 
never very small in the present approach (in contrast 
r ~ .11 to .042 in PM theory). Hence, the ferromag- 
netic exchange may be treated within Stoner mean-field 
theory, with no significant effect of the long-wavelength 
PM fluctuations. This justifies a posteriori our treat- 
ment of this coupling within the large-connectivity limit. 
At low-energy, we have a liquid of quasi-particles char- 
acterized by the effective Fermi scale e* F . In the ab- 
sence of any magnetic exchange (I = 0), the suscepti- 
bility of this gas would be of order Xqv — V e F- The 
actual susceptibility x — 1/( £ f — I) is correctly given 
by Stoner expression, with an effective Stoner enhance- 
ment S e ff = x/Xqp — -G6Rw- S e ff depends weakly on 
pressure, and measures the fraction of the total suscep- 
tibility enhancement due to the exchange (in contrast, 
S = x/xa = l-5S e f f(m* /m) is a combination of ex- 
change and localization effects and strongly depends on 
pressure). These remarks also imply that there is no 
significant enhancement of the effective mass due to fer- 
romagnetic spin fluctuations (in contrast with the loga- 
rithmic effect of PM theory), which leaves the estimate 
m*/m ~ 1-1/(1 — U/U c ) used above essentially unaf- 
fected. From Eq.(|^) and the fact that I > 0, one sees 
that the susceptibility is peaked around q = 0. For low 
\q\ = q Uf and u> -C qv* F (with v* F — Zv F the effective 
Fermi velocity), we can approximate the dynamical sus- 
ceptibility by (neglecting all other residual interactions 
between quasiparticles apart from the exchange): 
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From this expression, we see that there is a spin- 
fluctuation peak in Imx at a frequency uj max (q) ~ qvF / S. 
The peak height is of order Imx maa; ~ S/cp. These esti- 
mates coincide with those found in conventional PM the- 
ory, and are in reasonable agreement with the available 
neutron data |14| ). In contrast, the correlation length of 
the ferromagnetic fluctuations is found in our approach 
to be of order IcfS, — y/Seff — \/Rw , and hence much 
shorter and less pressure dependent than in PM the- 
ory where kp^ ~ y/S. A direct study of this quan- 
tity would help to clarify the nature of spin fluctuations 
in liquid 3 He. We also observe that in the whole do- 
main where q -C kp, the q 2 /k F terms in EqJ| is always 
negligible in front of 1 — I /e* F . Hence, the spin fluc- 
tuation mode has always linear dispersion uj max oc q, 
in contrast to PM theory where a regime u> max oc q 3 
is reached above a characteristic wavevector oc kp/\fS 
and below a low-energy scale oc €f/S 3 ^ 2 . Here in con- 
trast, a single energy scale exists T F ~ e F ~ .GGRw^f / S 
(~ 700mK at p = 0, ~ 400mK at p = 34bars). As 
a result, the susceptibility is not predicted to display a 
1/T 4 / 3 dependence at high temperature as in PM the- 
ory but rather follows a Curie law \ hi better 
agreement with experiment. The low-temperature be- 
havior of the susceptibility however is, rather remark- 
ably, predicted to be quite similar to the PM result 
even though Stoner mean-field theory applies to the low- 
energy quasiparticles. Indeed, the latter yields a low- 
temperature correction: e* F x = S e ff[l - cS e ff{T/T F ) 2 ] 
so that: x/xo = S[l - c'/R w {ST/T F ) 2 }. This expres- 
sion obeys the scaling T X (T) = f{T/T* F ) (since R w (p) 
weakly depends on pressure) and is similar to the PM 
result which is in good agreement with experiments 
and differs from the result of naive Stoner theory. 

Finally, we briefly discuss the magnetic field depen- 
dence of the magnetization, which has been proposed 
as a way to discriminate between the "almost ferromag- 
netic" and "almost localized" approaches || . Recently, a 
remarkable experiment []l5| (at p — 26bars) has revealed 
that the metamagnetic behavior predicted within the GA 
H at h^ A ~ 25T is not observed, up to (effective) fields 
of the order of 200T. In the present model, the magne- 
tization as a function of field can be deduced from the 
corresponding result for the Hubbard model m = fji{h), 
by solving the equation: m = fu{h + Im). The func- 
tion fu has been computed numerically in a previous 
work Jl6{ . Using the parameters (U,I) corresponding 
to p — 26bars we find that model (|^) does present at 
half-filling a metamagnetic transition, though at a much 
higher field h c ~ 80T than the GA result. This indicates 
that the constraint of half-filling is too strict to provide 
a reasonable lattice-gas description of the magnetization 
experiments. We have checked that extending the model 



to allow for a small concentration of vacancies 8 ~ 8% 
(in the spirit of Ref. ||) allows a reasonable description 
of the experimental magnetization curve m(h). 

In conclusion , we have proposed a physical picture of 
normal 3 He as a "Mott-Stoner" liquid , which seems in 
qualitative agreement with several experimental apsects. 
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